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Nucleosynthesis of s-elements in zero-metal AGB stars 
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Abstract. Contrary to previous expectations, recent evolutionary models of zero-metallicity stars show that the 
development of mixing episodes at the beginning of the AGB phase allows low- and intermediate-mass stars 
to experience thermal pulses. If these stars, like their metal-rich counterparts, also experience partial mixing of 
protons from the H-rich envelope into the C-rich layers at the time of the third dredge-up, an extensive neutron 
capture nucleosynthesis leads to the production of s-process nuclei up to Pb and Bi. Nucleosynthesis calculations 
based on stellar AGB models are performed assuming a parameterized H-abundance profile below the convective 
envelope at the time of the third dredge-up. Despite the absence of Fe-group elements, the large neutron flux 
resulting from the ^^C (a , n) reaction leads to an efficient production of s-process elements starting from the 
neutron captures on the C-Ne isotopes. Provided partial mixing of protons takes place, it is shown that population 
III AGB stars should be enriched in s-process elements and overall in Pb and Bi. 
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1. Introduction 



The determination of chemical abundances at the sur- 
face of evolved stars provides important clues on their 
structural evolution and traces the hidden processes tak- 
ing place in their nuclearly active interior. Asymptotic gi- 
ant branch (AGB) stars, and in particular zero-metallicity 
stars, represent an ideal laboratory for such a study. They 
correspond to the late phase of the evolution of stars with 
masses between about 1 and 8 Mq , which includes more 
than 80% of all stars; they exhibit peculiar surface abun- 
dances as compared to other red giant stars; and many 
of them are characterized by a strong mass loss (up to 
10~^Mq yr~^) which eject the surface material into the 
interstellar medium, contributing thereby to the galactic 
chemical evolution. Among those, the population III stars, 
i.e the first to form with the chemical composition inher- 
ited from the Big Bang nucleosynthesis, are of particular 
interest to understand the early enrichment of the uni- 
verse in elements heavier than hydrogen and helium. Little 
is known observationally on population III stars, since 
so far no star of population III initial composition has 
yet been found. If s-process enriched AGB stars are com- 
monly observed nowadays, they all correspond to metal- 
licities [Fe/H] ^ —3. However, the detection of carbon 
and nitrogen enrichments in the intergalactic medium at 
high redshifts and in a considerable number of extremely 
metal-poor stars of our Galaxy suggests the existence of 
primordial AGB stars (see e.g Abia et al. 2001). 

The abundance peculiarities observed at the surface 
of AGB stars are understood as resulting from the mix- 



ing of material synthesized in the interior with the surface 
layers as the convective envelope deepens during the so- 
called third dredge- up (denoted hereafter 3DUP). But the 
nucleosynthesis of s-process elements requires neutron-rich 
environments that cannot be accounted for by the sole ac- 
tion of the dredge-up. One way to solve the problem is 
to assume the partial mixing (PM) of protons from the 
envelope into the ^^C-rich layers during a 3DUP event 
(Iben & Renzini 1982, Busso et al. 1999). The incomplete 
operation of the CN cycle forms a ""'^^C pocket" which 
burns radiatively and releases a large number of neutrons. 
However, AGB models are still subject to large uncertain- 
ties concerning the consistent prediction of both the 3DUP 
and PM processes (e.g Mowlavi 1999). In particular, the 
3DUP and PM properties are sensitive to the stellar char- 
acteristics (such as the mass, metallicity or mass loss rate) 
and to the numerical scheme used to describe the convec- 
tive boundaries (e.g Frost & Lattanzio 1996). In the past 
years, the mixing scenario has been revived by the incor- 
poration of new physical ingredients to model overshoot- 
ing (Herwig et al. 1999) and rotation (Langer et al. 1999). 
Models including these new treatments predict the forma- 
tion of a ^•^C pocket in a zone comprising some 1O~^M0. 
However, the total amount of ^^C produced in this way de- 
pends on the PM parameters (extent and efficiency of the 
overshoot, amplitude of rotation) which cannot be quan- 
tified ab initio at the present time. In this context, the 
resulting uncertainties on the subsequent nucleosynthesis 
can only be handled at the moment in a parametric way. 
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Recently several works have been devoted to the study 
of the structure and evolution of primordial stars (Siess et 
al. 2001, Chiefh et al. 2001, Marigo et al. 2001). However, 
from a nucleosynthesis point of view, the production of el- 
ements heavier than Mg and ultimately of iron in low- and 
intermediate-mass population III stars remains largely a 
mystery. Assuming that primordial AGB stars suffer PM 
of protons as their higher metallicity counterparts do, we 
analyze its effect on the resulting nucleosynthesis applying 
the PM model developed by Goriely and Mowlavi (2000) 
to the realistic 1-D AGB models computed by Siess et al. 
(2001). The corresponding modeling is described in Sect. 
2. In Sect. 3, we study in detail the nucleosynthesis oper- 
ating in the PM zone and its consequences on the surface 
abundance predictions of zero- metallicity AGB stars. In 
particular, it is shown that the initial absence of Fe does 
not prevent the formation of the highest-mass elements 
and leads to the overproduction of Pb and Bi. 

2. Astrophysical model and input physics 

The stellar models used in this study arc part of a re- 
cent grid computed by Siess et al. (2001) and we refer 
the reader to this publication for ampler details. Among 
the available mass tracks, we selected the 3 Mq model for 
two main reasons: it is our most advanced model which 
accounts for 22 thermal pulses and 3M0 is a somewhat 
"representative" mass for an AGB star, between the low- 
(IM0) and massive (8M0) ones. Briefly, the evolution 
of this model (which lasts for about 250 Myr) proceeds 
as follows: during the main sequence, the star cannot 
activate the CNO cycle because of the absence of ^^C. 
Consequently, it contracts and develops higher central 
temperatures to sustain its energy budget by the poorly 
efficient p-p chains. When the central temperature reaches 
lO^K, the 3a reactions start and carbon is produced in the 
core. The CNO cycle finally activates when the ^^C mass 
fraction exceeds ~ 10~^^. At the exhaustion of hydrogen 
in the center, the 3a reactions are still effective and the 
star enters the core He burning phase without experienc- 
ing the first dredge-up. The beginning of the AGB phase 
of zero-metallicity stars is characterized by the occurrence 
of mixing episodes which enables the carbon produced in 
the He burning shell to mix with the envelope. This car- 
bon pollution of the structure activates the CNO cycle in 
the H burning shell and allows primordial stars to resume 
a "standard" AGB evolution comparable to that of more 
metal rich stars. Note that without this mixing episode, 
the 3 Mq models would probably not suffer thermal pulses 
(and 3DUP). Primordial AGB stars are characterized by 
massive CO cores as a result of initially expanding their 
nuclear core to compensate for the non- activation of the 
CNO cycle. They also rapidly undergo third dredge-up 
events (after 9 pulses in this case) and exhibit higher tem- 
peratures in the pulse compared to more metal-rich stars. 
Of nucleosynthetic relevance is the activation of very ef- 
ficient hot bottom burning (HBB) after the 20th pulse, 
where the temperature at the base of the convective en- 



velope reaches ~ 8 x 10''K. The impact of HBB on the 
surface abundances of CNO elements is discussed in de- 
tails in Siess et al. (2001). Finally, because of the low 
metal content in the envelope, mass loss is rather ineffi- 
cient and by the end of the computations only 1O~^M0 is 
lost. Nevertheless, the successive dredge up episodes will 
increase little by little the opacity in the envelope and 
consequently the mass loss efficiency. 

The temperature and density profiles given by the stel- 
lar evolution code are used as input data to compute the 
detailed nucleosynthesis resulting from the PM of protons. 
Protons from the envelope are ingested artificially from 
the bottom of the convective envelope into the underlying 
C-rich region at the time of the 3DUP events. We adopt 
two different H-abundance profiles. The first one decreases 
exponentially with depth and ranges from X™^^ = 0.7 
at the bottom of the convective envelope to 10~^ at the 
bottom of the PM zone. The second H abundance pro- 
file decreases slower with depth than the previous case. It 
corresponds to the "Slow" curve in Fig. 10 of Goriely and 
Mowlavi (2000). However, since the nucleosynthesis calcu- 
lations are found to be almost identical in both cases, we 
present here only results using the second H abundance 
profile which can be considered as more characteristic of 
diffusion processes (Herwig et al. 1999). We consider a PM 
zone of extension Apm = 0.05 corresponding to 5% in mass 
of the pulse-driven convection zone. This corresponds to 
a zone comprising some 10~^ Mq subject to the s-process 
neutron irradiation. 

The nucleosynthesis resulting from the ingestion of the 
protons into the C-rich layers at the time of each 3DUP is 
followed during the first 22 interpulse and pulse phases 
of the 3Mq stellar model. In contrast to Goriely and 
Mowlavi (2000), the ashes of the interpulse and pulse nu- 
cleosynthesis are reingested in the next interpulse/pulse 
sequence, so that multiple irradiations on former nucle- 
osynthetic products are consistently taken into account. 
The nuclear reaction network includes 817 nuclei up to Po, 
and all relevant nuclear, weak and electromagnetic inter- 
actions. Note that the network needs to be extended to the 
neutron-rich region because of the very large neutron den- 
sities generated in the late He-flashes. The latest experi- 
mental and theoretical reaction rates are taken from the 
Nuclear Astrophysics Library of the Brussels University 
(available at http://www-astro.ulb.ac.be). More details 
can also be found in Goriely and Mowlavi (2000). 



3. The neutron capture nucleosynthesis 

To analyse the nucleosynthetic impact of the PM of pro- 
tons into the C-rich region, we first consider one inter- 
pulse/pulse/3DUP phase. The richest nucleosynthesis oc- 
curs during the interpulse phase. More specifically, for low 
proton-to-carbon ratios [Yp/Y{^^C) ^ 1, where Y{A) = 
X{A)/A is the abundance of A by number], neutrons are 
efficiently produced after -'^^C has been synthesized, and 
subsequently destroyed by ^'^C (a , n) ^^O (Fig. |l|) . In the 
absence of Fe-group nuclei, neutrons are captured by all 
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the abundant C, N, O, F and Ne isotopes and enable the 
production of high mass elements. The fast (n, a) reac- 
tion on ■^'^S forms a major bottle-neck braking down the s- 
process flow partially. As soon as the bottle-neck is passed, 
radiative neutron captures followed by /3-decays on unsta- 
ble nuclei lead to the synthesis of all medium and heavy 
elements up to Pb and Bi. The ^^S (n , a) ^°Si recycling 
reduces the total amount of s-elements produced by a fac- 
tor of about 5. The evolution of the mass fraction of some 
representative nuclei (with a magic number of neutrons 
TV = 28, 50, 82 and 126) is displayed in Fig. |l|. In the se- 
ries of neutron captures and /3-decays characteristic of the 
s-process, iron group nuclei do not play any particular 
role, in contrast to the classical picture of the s-process 
where the neutron captures take place on the abundant 
^^Fe mainly. 
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Fig. 1. Mass fraction of some elements representative of 
the neutron capture nucleosynthesis taking place in the 
PM layer characterized by an initial proton-to-^^C ratio 
of 10~^ during the 21st interpulse phase of the M — 
3 Mq,Z = AGB star. Also given (bold sohd line) is 
the time dependence of the neutron exposure t. 

For 10^^ ^ Yp/Y{^^C) ^ 1, the neutron exposure 
T — Jq Nn VT dt reaches 10 mbarn"^ {vt is the most 
probable relative neutron-nucleus velocity at tempera- 
ture T) and the neutron density A^n a few 10^ cm~^. 
These high values explain the large production of met- 
als up to Pb and Bi despite the absence of initial iron 
(Fig. Ija). For relatively high proton-to-^^C abundance ra- 
tios [Yp/Y{^^C) ^ 1], i-^N is more abundant than ^^C 
and prevents any significant neutron production. In this 
case, the PM of protons only affects light species. The 
fast ^''N (n , p) ^^C reaction is responsible for the pro- 
duction of ^^C and later of by {a , 7) ^^O . This 
reaction also releases a large amount of protons which 
are captured by ^^O to form ^^N. Subsequently, in the 
pulse-driven convection zone, ^^F is largely produced by 
^^N {a , 7) ^^F . However, as soon as envelope HBB is acti- 
vated, the dredged-up ^^F is destroyed. HBB also increases 
the surface abundance of ^''N and ^^Na significantly at the 
expense of CNO isotopes and ^^Ne, respectively. 
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Fig. 2. a) Abundance distribution (relative to solar) at the 
end of the 21st interpulse phase (a) and following the next 
pulse (b) for the M = 3 M©, Z = AGB star. In panel (a), 
the distribution is averaged over the mass range covered 
by the PM region (in this case Mpm = 4 x 10"** Mq) while 
case (b) gives the surface abundance after the 3DUP. Full 
squares correspond to s-only isotopes. Note the nuclear 
rearrangement of some s-only elements during the pulse. 

As the core mass grows, the temperature in the ther- 
mal pulse increases and by the 20th pulse, it exceeds 
4 X lO^K. In these conditions, ^^Ne is rapidly burnt by 
^^Ne {a , n) ^^Mg and additional neutrons are released. 
In the 22nd pulse, the neutron density reaches A^n = 
10^^ cm~'^, but for a short period of time and mostly 
near the bottom of the He-flash zone. Consequently, 
the pulse-averaged resulting neutron exposure is small. 
Nevertheless, the neutron burst allows for some rearrange- 
ments of the abundance distribution of peculiar s-only 
isotopes (e.g ®°Kr or ^"^^Nd, see Fig. ^) and of some sr- 
nuclei by short-cutting the s-process flow. In particular, 
^^^Sn, ^'^^Xe, ^^^Ba and ^"^^Ce are overproduced compared 
with neighboring isotopes. Abundances of the so-called T- 
dependent s-branching nuclei, e.g '^^'^Er, are also modified 
because of the high temperatures prevailing in the pulse. 
Among the medium-mass nuclei, ^^S, ''"Ar and "^^Ca are 
largely produced. Note that the large dip in the X/Xq 
scale observed in Figs. |§-||for the Fe-group elements is due 
to their high solar values, and consequently small abun- 
dances relative to solar. 

Except for the first 6 interpulse/pulse sequences where 
the C abundance left over by the preceding pulse-driven 
convective zone is small, the nucleosynthesis resulting 
from the PM of protons is found to be qualitatively iden- 
tical from one interpulse/pulse phase to the other (the 
only significant difference originates from the tempera- 
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tures reached at the bottom of the He-flash zone). The 
relatively small extent of the PM zone (Apm — 0.05) does 
not allow a significant cumulative neutron irradiation on 
the ashes of the previously synthesized material. As a con- 
sequence, the most abundant material, namely the ele- 
ments comprised between C and Ne, always remain the 
seed nuclei for the neutron captures in the PM zone. 

The above-described nucleosynthesis is in many re- 
spects similar to the one found in the low-metallicity 
{Z = 0.001) star subject to the same mixing of protons 
(Goriely & Mowlavi, 2000). In the radiative C-rich layers 
during the interpulse phase, the overall abundance distri- 
butions are much the same even if the seed nuclei at the 
origin of the heavy elements are in both cases quite dif- 
ferent. For stars with metallicities Z ^ 0.001, the number 
of neutrons available per seed nuclei exceeds the number 
needed to produce the heaviest s-nuclei, i.e Pb and Bi, so 
that the type of seeds as well as their initial abundances 
only affect the absolute abundance of heavy elements syn- 
thesized by the s-process, but not their relative distri- 
bution. However, for extremely low metallicity stars, the 
higher temperatures reached at the bottom of the pulse- 
driven convective zones can reshape the s-abundance dis- 
tribution differently. 

Provided partial mixing of protons occurs in popu- 
lation HI stars in a similar way as in more metal-rich 
stars, the primordial stars should exhibit a strong over- 
abundance of Pb and Bi compared with other s-process 
elements ([Pb/s] ^ 1). Therefore, such stars could be 
dubbed population HI Pb-stars (Goriely & Mowlavi 2000, 
Van Eck et al. 2001). In addition to their large overabun- 
dance in ^^N (mainly due to the HBB transmuting all 
the CNO envelope isotopes into ^''N), these Pb-stars are 
also predicted to have a large ^'^Na content, as shown in 
Fig. ^. Van Eck ct al. (2001) observed low-metallicity stars 
(down to [Fe/H]=-2.45) enriched in s-process elements and 
characterized by large Pb overabundances ([Pb/hs]J; 1, 
where hs denotes heavy s-elements such as Ba, La or Ce). 
The discovery of these Pb-stars confirms the s-process 
efficiency in low-metallicity stars. However, this conclu- 
sion is not reached by Aoki et al. (2000) who found in a 
metal-poor s-process enriched stars ([Fe/H]=-2.7), a Pb 
overabundance ([Pb/Fe]=+2.6) but no Pb-to-s overabun- 
dance ([Pb/hs]~0). This latter observation remains to be 
explained. 

While the discovery of Pb-stars (Van Eck et al. 2001) 
strongly confirms the "proton-mixing" scenario for the de- 
tailed operation of the s-process in AGB stars, their is still 
no evidence neither on the existence of population HI stars 
nor on the possible mixing of protons in such stars. For 
this reason, it would be premature at this stage to esti- 
mate the possible contribution of population HI stars to 
the galactic enrichment in s-process elements, and more 
particularly in Pb (e.g Travaglio et al. 1999). 
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Fig. 3. Elemental abundances (full circles) and overabun- 
dances relative to solar (open circles) at the surface of 
the M — 3 Mq zero-metallicity AGB star after the 22nd 
convective pulse phase. 

4. Conclusion 

Assuming that population HI stars experience, just like 
more metal-rich stars, a partial mixing of protons from 
the H-rich envelope into the C-rich layers at the time 
of the third dredge-up, an extensive neutron capture nu- 
cleosynthesis is found to take place leading to the effi- 
cient production of s-process nuclei. Despite the absence 
of Fe-group elements, the large neutron flux resulting from 
the ^^C (a , n) ^^O reaction transmutes part of the C-Ne 
isotopes into s-nuclei with a large overabundance of Pb 
and Bi. Provided primordial stars are observable today, if 
they exhibit a surface enrichment in s-process elements, a 
strong overabundance of Pb and Bi compared with other 
s-process elements ([Pb/s] ^ 1) is predicted. These pop- 
ulation HI Pb-stars remain to be observed. Their discov- 
ery or the indirect identification of their pollution in later 
stellar generations would prove that the galactic enrich- 
ment in s-process elements could have started already a 
few hundred million years after the Big Bang. 
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